Absence of the d-Density Wave State in 2D Hubbard Model 
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Using the Dynamical Cluster Approximation (DCA) we calculate the alternating circulating- 
current susceptibility and investigate the transition to the d-density wave (DDW) order in the 
two-dimensional Hubbard model. The 2x2 cluster used in the DCA calculation is the smallest that 
can capture d-wave order; therefore, due to the mean-field character of our calculation, we expect 
to overestimate d-wave transition temperatures. Despite this, we found no transition to the DDW 
state. In the pseudogap region the DDW susceptibility is enhanced, as predicted by the slave boson 
SU(2) theory, but it still is much smaller than the d-wave pairing susceptibility. 
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Introduction The high T c cuprates display a vari- 
ety of unusual properties, which remain unexplained by 
conventional theories. The most intriguing physics oc- 
curs at small doping, in the proximity of antiferromag- 
netism and superconductivity, and is characterized by 
non-conventional behavior of many observables includ- 
ing the spin susceptibility, optical conductivity, specific 
heat and transport properties. Many of these unusual 
properties are associated with the presence of a pseudo- 
gap in the one- and two-particle spectra. Photoemission 
spectra show that, at small doping and above T c , the 
states around the (tt, 0) point in the Brillouin Zone are 
gapped and Fermi segments appear around (7r/2,7r/2), 
suggesting that the symmetry of the pseudogap in the 
hole-doped cuprates is consistent with the d-wave sym- 
metry of the superconducting gap 0. 

Based on the cuprate phenomenology, Chakravarty et 
al. 0] proposed that the pseudogap results from the com- 
petition between two ordering processes. One is d-wave 
superconductivity (DSC) and the other is a state char- 
acterized by long-range order of alternating orbital cur- 
rents. The latter is a staggered-flux state which breaks 
the translational and the time-reversal symmetry and 
represents in fact a charge density with d-wave symme- 
try Eli i- e - it is a d-density wave (DDW) state. In this 
scenario the system evolves continuously from the DDW 
state to the DSC state with decreasing temperature or 
increasing doping, and the two states coexist up to op- 
timal doping. The experimentally observed one-particle 
spectra in the pseudogap region can be well understood 
on the basis of the DDW state pj], which makes it a 
very appealing candidate for the origin of the pseudogap 
physics. Other properties of the cuprates, such as the 
resonant peak in the superconducting state and the dop- 
ing dependence of superfluid density seem also to be well 
captured by this model Recently it was proposed the 
interplay of the DDW and the inter-planar tunneling of 
Cooper pairs to be responsible for the T c dependence on 
the number of CuC>2 layers which characterizes different 
materials p- In principle, the presence of the DDW 
state has subtle experimental consequences, such as the 
formation of a magnetic moment associated with the or- 



bital currents Q • The interpretation of the experimental 
data in this respect however is still controversial. 

Whereas the theory of the DDW is phenomenological, 
slave boson theory holds the promise of a microscopic 
basis which may explicitly consider static or fluctuating 
DDW order, as well as d-wave pairing. These are uncon- 
trolled theories for the t — J model, which is equivalent 
to the strong coupling limit of the Hubbard Hamilto- 
nian. They are influenced by P.W. Anderson's et al. 
idea of resonance valence bond state. The charge and 
the spin degrees of freedom are separated by introducing 
auxiliary slave bosons. The resulting mean- field theories 
explicitly decouple the fermion hopping along the bonds, 
the fermion pairing and the bosonic field, and produce 
phase diagrams similar to the experimental one. 

We briefly discuss the main results of slave boson theo- 
ries. The t — J Hamiltonian has a local SU (2) symmetry 
at half filling . As a result of this symmetry the 7r-flux 
state (a staggered- flux state with the flux per plaque- 
tte equal to tt) and the d-wave paring state are degener- 
ate in the undoped model. Doping breaks the symmetry 
to U(l) and the d-wave state becomes energetically fa- 
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vored [lfj. The d-wave state is characterized by a finite 
value of the fermion pairing operator and a real super- 
conductor emerges below the condensation temperature 
of the bosons. At low doping (S < 0.05) and for T > 
the d-wave pairing becomes unstable towards the 7r-flux 
or staggered-flux state of spinons [n). This is the stan- 
dard picture of the Z7(l) slave boson mean field theory of 
the t — J model. However, the inclusion of other terms in 
the mean- field decoupling, such as the holon's flux [l2j |. 
results in a existence of a DDW state at finite doping and 
above T c , but exclude the coexistence of the two states. 

One of the drawbacks of the U(l) theory is that its so- 
lution is not stable against the fluctuations of the gauge 
field H3- Fluctuations are especially important at small 
doping, where the energy difference between states con- 
nected via a SU (2) transformation is very small (since 
they are degenerate at zero doping). Therefore all these 
states have an important contribution in the determina- 
tion of the free energy. 
P.A. Lee et al. developed a slave boson mean field 
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theory which is SU{2) symmetric at finite doping |l4j . 
The price paid is that one must deal with two slave 
bosonic fields and three constraints. The advantage of 
this approach is that the SU{2) mean-field solution is 
likely superior at small doping, since it accounts better 
for the fluctuations between different low energy SU{2) 
connected states. Their solution for the pseudogap re- 
gion is a staggered- flux of fermions which is gauge equiv- 
alent with the ei-wave pairing of fermions. However, 
the fcrmion staggered-flux state is not the same as the 
staggered- flux state of electrons (or the DDW state), 
and neither breaks time-reversal nor translation symme- 
try. Therefore, in SU(2) theory, the pseudogap is not 
a broken symmetry state with long range order as it is 
DDW, but is rather is characterized by strong spatial and 
dynamic fluctuations between d-wave, s-flux and other 
SU(2) related states. 

The goal of this paper is to investigate the inter- 
play between DDW and DSC order in the 2D Hubbard 
model. U sing the Dynamical Cluster Approximation 
(DCA)[H Eg we calculate the response functions as- 
sociated with these two types of order. The DCA sys- 
tematically adds non-local corrections to the Dynamical 
Mean Field Approximation (DMFA) [ril flil ] by mapping 
the lattice onto a finite-size periodic cluster. The DCA 
mapping from the lattice to the cluster is accomplished 
by coarse-graining all the internal propagators in irre- 
ducible Feynman graphs in reciprocal space. Correlations 
at short length scales, within the cluster, are treated ex- 
plicitly with a quantum Monte Carlo (QMC) simulation, 
while those at longer length scale are treated at the mean 
field level. Due to the residual mean field character of 
our approximation, we expect our calculation to overes- 
timate the transition temperatures of both the DDW and 
DSC critical temperatures. 

Generally, we expect to see the the most pronounced 
mean-field behavior from the smallest cluster that can 
reflect the broken symmetry. A similar situation occurs 
in DMFA simulations of the antiferromagnetic phase of 
the Hubbard model, where Neel order is possible since 
the impurity spin and the mean-field host may have op- 
posite spin orientations. Since non-local fluctuations are 
suppressed, the DMFA overestimates the Neel transition 
temperature. In the present case, the N c — 4 (i.e. 2x2) is 
the smallest possible cluster allowing for <i-wave pairing 
or a circulating current. Orbital antiferromagnetism is 
possible since the moment in the cluster and the host can 
have opposite orientations. Since fluctuations on longer 
length scales are suppressed, we would expect to overesti- 
mate both the d-wave superconducting and the d-density 
wave transition temperatures. 

Formalism We present DCA calculations for the con- 
ventional 2D Hubbard model describing the dynamics of 
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FIG. 1: The one particle total and K dependent DOS at 
S = 0.05 doping. Inset: The uniform magnetic susceptibility 
versus T. The maximum defines the pseudogap temperature 
T*. 



electrons on a square lattice. The model 
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is characterized by a hopping integral t between nearest 
neighbor sites and a Coulomb repulsion U two electrons 
feel when residing on the same site. As the energy scale 
we set t = 0.25eV so that the band-width W = 8t = 2eV, 
and study the intermediate coupling regime U = W . 
We study the dynamics on short length-scales by set- 
ting the cluster size to N c — 4. This cluster size is large 
enough to capture the qualitative low-energy physics of 
the cuprate superconductors [l^. Eof . The correspond- 
ing phase diagram resembles the generic phase diagram 
of cuprates |20( , displaying regions characterized by anti- 
ferromagnetism, (i-wave superconductivity, Fermi liquid 
and pseudogap regimes, in qualitative agreement with 
experimental results. 

In this paper we calculate the static (i.e. u> = 0) sus- 
ceptibilities which correspond to the circulating current 
(cc) operator, 



W 



(2) 



and respectively to the d-wave pairing operator, 

P = E-9( fc ) Cfe l c - fc T > ( 3 ) 

k.a 

where g(k) — cos{k x ) — cos(k y ) is the d-wave symmetry 
factor. 

Results The pseudogap temperature, T* is deter- 
mined from the maximum in the uniform magnetic sus- 
ceptibility (see the inset in Fig. ^) when accompanied by 
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FIG. 2: The d-wave pairing (circles) and the circulating cur- 
rent (squares) susceptibilities versus temperature, at a) 5% 
and b) 25% doping. The dotted vertical line in a) is at the 
pseudogap temperature T* . 

a suppression of spectral weight in the DOS. We show 
this in Fig.^where the total and the fC-dependent DOS, 
below T*, at 8 = 0.05 doping is plotted. The DC A on a 
N c = 4 cluster implies a coarse graining of the Brillouin 
Zone in four cells around K = (0,0), (0, ir), (tt,0) and 
(tt, 7r) and the if-dependent DOS corresponds to the av- 
erage over all k belonging to a coarse-grained cell of the 
single particle spectra A(k,u>). This poor resolution in 
the reciprocal space allows to study only the gross fea- 
tures in the single-particle spectra. Despite this, it can 
be seen from Fig.^that the pseudogap in the total DOS 
is a result of the suppression of spectral weight in the 
cell at (0, tt). Therefore, we believe that our calculations 
capture well the experimentally observed features of the 
pseudogap. 

However, our calculations show that these features are 
not a consequence of the DDW state. In Fig. |2J-a we 
plot both the cZ-wave pairing susceptibility and the cc- 
susceptibility versus temperature, at S — 0.05 doping. 
The pairing susceptibility diverges at T c , indicating a d- 
wave superconducting instability. The cc-susceptibility 
does not diverge, indicating the absence of a possible 
transition to the DDW state. 

At large temperatures the d-wave pairing and the cc- 
susceptibilities are degenerate, and they both increase 
with decreasing temperature. In the pseudogap region 
(left side of dotted line) the eZ-wave pair field susceptibil- 
ity is much larger than the cc-susceptibility. Close to T c , 
the cc-susceptibility saturates and starts even decreasing 
with decreasing temperature. The fact that in the pseu- 
dogap region both cZ-wave and the cc-susceptibilities are 
enhanced show that fluctuations between these states are 
significant, as it was predicted by the SU(2) theory [lij . 

In Fig. |21-b we show the d-wave and the cc- 
susceptibilities at S = 0.25 doping. No pseudogap [2(J 



is observed at this doping. We notice that, starting well 
above T c , the cc-susceptibility decreases with decreas- 
ing temperature. This behavior is different from the one 
observed at small doping where the cc-susceptibility in- 
creases with decreasing T up to T c . We therefore con- 
clude that in the overdoped region the fluctuations be- 
tween DSC and DDW above T c are much less important. 

Calculations (not shown here) with other values of the 
parameters (different values of U), or with the inclusion 
of next-nearest-neighbor hopping corresponding to both 
electron and hole doping, exhibit similar results. A diver- 
gent DDW susceptibility is never found. Our results are 
consistent with renormalization group studies |2l| where 
no divergence of the cc-susceptibility is found, and with a 
mean-field treatment of an extended Hubbard mo del [22]] 
where an additional correlated hopping term was neces- 
sary to stabilize the DDW state. 

Conclusions We present a DCA calculation of the 
two-dimensional Hubbard model, focusing on the compe- 
tition between DDW and DSC orders. We showed pre- 
viously that the DCA calculation captures the generic 
features of the pseudogap region as seen in the photoc- 
mission and magnetic measurements. Nevertheless, as we 
show here, these properties are not a consequence of the 
existence of DDW state. The DCA should overestimate 
any DDW transition temperature but, despite this, we 
found no transition to such state. 

We also found that both the cc-susceptibility and the 
eZ-wave pairing susceptibility are enhanced in the pseu- 
dogap region, indicating that the fluctuations between 
these states is significant. This is not true in the over- 
doped region, where we found that the cc-susceptibility 
is suppressed above T c . 
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